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VAPORIZATION RATES AND DRAG COEFFICIENTS FOR ISOOCTANE 
SPRAYS IN TURBULENT AIR STREAMS 


By Robert D. Ingebo 


SUMMARY 


A droplet camera developed at the NACA Lewis laboratory was used 
to obtain drop-size distribution and drop-velocity data for isooctane 
injected from a simple orifice directly into a turbulent air stream. 
From these data and wet-bulb temperature data, vaporization rates and 
drag coefficients were calculated for isooctane drops accelerating and 
evaporating in streams having velocities of 140 and 180 feet per second. 
Spray vaporization rates based on the mean drop diameter Dog were com- 


pared with single-drop vaporization rates by use of the heat-transfer 
equation 


ар2 o Ak At 
E i E. А 
99 рун, 


where 020 is defined ав the diameter of a drop having an area equal to 


the ratio of total area to total number of drops formed in the spray, Ө 


is the time, kg is the thermal conductivity of the air, At is the dif- 


ference between the air temperature and drop-surface (wet-bulb) temper- 
ature, рт, ів the density of the liquid, Н. is the latent heat of 


vaporization, and Nusg is the heat-transfer Nusselt number based on 
the mean drop diameter Dog. For these tests, the change in sensible 


heat of the liquid was negligible compared with latent heat requirements. 
The weight percent of the spray evaporated at distances of 1 to 18 in- 
ches downstream of the injector was determined and found to be in agree- 
ment with data obtained in a previous investigation by means of a 
sampling-probe technique. 


Nukiyama-Tanasawa and log-probability drop-size-distribution func- 
tions gave good agreement with experimental drop-size-distribution data 
in the determination of mean drop size. The Sauter mean diameter for the 
spray, predicted by the Nukiyama-Tanasawa empirical expression for a 
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different type of nozzle, agreed with that found for the 180-feet-per- 
second air stream. The maximum drop diameter for the isooctane sprays 
was calculated from a balance of aerodynamic pressure and surface ten- 
sion forces under the condition of maximum acceleration of the sprays. 
These values compared well with the maximum drop diameter observed from 


photomicrographs of the sprays. 


INTRODUCTION 


When а: liquid fuel is injected into the preignition zone of а ram 
jet or a turbojet engine, or into an afterburner, the liquid is broken 
up into a cloud of droplets which are then accelerated to stream veloc- 
ity. While accelerating, they evaporate at a rate determined by: (1) 
the air-stream velocity, temperature, and static pressure, (2) the 
droplet velocity, temperature, and diameter, and (3) the physical prop- 
erties of the liquid and vapor. Thus, the vaporization rate of sprays 
in air streams of known temperature, pressure, and velocity may be de- 
termined by use of heat-transfer equations when data on drop-size dis- 
tribution, drop acceleration, and drop-surface temperature are available. 


Several investigators (e.g., refs. 1 and 2) have obtained drop- 
size-distribution data for the atomization of liquids in air streams. 
However, the problem of correlating drop-size-distribution parameters 
with physical characteristics of the atomization process is still un- 
solved. Data have been obtained (refs. 3 and 4) for the weight percent 
of fuel spray evaporated in air streams at given distances downstream 
of the injector by means of & sampling technique. Since this method 
did not give data on drop size and drop velocity, a correlation with 
characteristics of momentum, mass, and heat transfer could not be made. 
Correlations of this type have been made for single drops evaporating 
in air streams (refs. 5 and 6). However, the problem of relating 
vaporization rates of single drops to that of sprays has not been solved 
because of the lack of spray vaporization-rate data. 


In reference 7, theoretical equations are given for predicting the 
vaporization rate of fuel spray based on a single mean-diameter droplet. 
For simplicity, an evaporation constant, assumed to be independent of 
drop size and velocity, was used. An expression valid only for either 
exceedingly fine sprays or sprays which have accelerated to stream ve- 
locity was obtained. If the doubtful assumption is made that evapora- 
tion prior to obtaining either of these conditions is negligible, the 
equation of reference 7 may be applied directly to the evaporation of 
fuel sprays in jet engines. | 


Та this investigation, data on drop-size distribution, drop veloc- 
ity, and drop-surface (wet-bulb) temperature were obtained for isooctane 
sprays in turbulent air streams having velocities of 140 and 180 feet 
per second. This was accomplished by using a droplet camera developed 
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at the NACA Lewis laboratory which gave photomicrographs of the sprays 
and also drop-velocity data.  Drop-size-distribution data were ob- 
tained from photomicrographs of the spray; and the wet-bulb temperature 
for isooctane, under the conditions of this investigation, was obtained 
by the experimental technique described in references 5 and 6. From 
these data, it was possible to calculate the following spray parameters: 
(1) mean-drop sizes,.(2) vaporization rates, (3) drag coefficients, and 
(4) weight percent evaporated. 


SYMBOLS 


The following symbols are used in this report: 


A surface area of drop, sq cm 
А! frontal area of drop, sq cm 
а, acceleration of drop, ft/ Sec? 
b - constant 
bg дг molecular mass diffusivity, g/(cm)(sec) 
Ср drag coefficient 
с root-mean-square velocity of air molecules, cm/sec 
D drop diameter, cm 
AD drop-diameter size range, cm 
D constant 
ЖарУ 
Doo; 2850» mean drop diameter defined by the expression рак = ape’ 
Dais- 32 for example, Doo is the diameter of a single drop hav- 


ing an area equal to the total area to number ratio for 
all the drops formed in а spray, cm 


E portion of total spray which has evaporated, percent 
F total force of air stream, dynes 
G momentum 


g acceleration due to gravity, 980 сш/ вес? 


п ~ P 


Ке 


At 
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latent heat of vaporization, g-cal/g 
heat-transfer coefficient, g-cal/(sec)(sq сш) (°C) 
thermal conductivity, g-cal/(sec) (sq em) (°c) /cm) 
distance from mirror to film, ft 

mean free path of air molecules, cm 
magnification factor for lens system, Zl power 
mass of liquid, g | 

vaporization rate, g/sec 

mass of liquid vaporized, g 

heat-transfer Nusselt number, hD/k 


heat-transfer Nusselt number based on mean drop diameter 
Doo» hDoo/k 


number of drops in given size range 


air-stream static pressure, in. Hg abs 
fuel injection pressure drop, 1b/sq in. 


volumetric flow rate, сс/вес 

constant 

volume fraction of drops having diameter >D 
Reynolds number based on drop diameter, D^up/u 
temperature, ОЕ 

temperature, 9C 


difference between air temperature and surface temperature 
of drop, t4 - tg, € 


velocity, ft / sec 
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Au 


We 


С 
б) 


Subscripts: 


relative velocity of air with respect to drop, cm/sec 
volume fraction of drops having diameter < D 


critical Weber number, p,U4D,,/20 


distance downstream from the injector orifice, in. 
in(D/D) 

constant 

vaporization time, sec 


evaporation constant, 1k. At/p E, sq сш/вес 


fluid viscosity, poises 
fluid density, g/ cu cm 
surface tension, dynes/cm 


mirror speed, rpm 


air 

drop 

final condition 
gas 

initial condition 
liquid 

maximum 

air stream 


vapor 
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APPARATUS AND PROCEDURE 


The vaporization rates of isooctane sprays evaporating and acceler- 
ating in turbulent air streams were determined experimentally with the 
apparatus shown in figure l. Air at 820559 Е and 1641 percent relative 
humidity was supplied from the central laboratory system at 40 pounds 
pe» square inch gage. ‘The air was metered with а variable-area orifice 
and exh&ugz! d to the altitude exhaust system. In order to minimize 
approach-stream turbulence, а 30 by 28 mesh, O.013-inch-diameter monel 
wire screen was placed 5 feet 4 inches upstream of the inlet to the 
transparent 8-inch-inside-diameter plastic test section. Two trans- 
parent 2- by 2- by 3/8-inch optically flat plastic windows were installed 
ir she test section for photographing the spray.  Air-stream pressure 
and velocity in the test section were controlled by upstream and down- 
stream valves. 


Isooctane fuel (2, 2, 4-trimethylpentane, A.S.T.M. specifications) 
was injected by pressurized nitrogen into the test section through & 
1/4-inch-diameter Inconel tube sealed at the end. А 0.041-inch-diameter 
orifice was drilled and reamed to size 1 inch from the sealed end in 
order to obtain a smooth surface. The orifice was pointed directly into 
the air stream and positioned on the test-section center line (position 
4, fig. 1) and also moved vertically to positions 1, 2, 6, and 7 in fig- 
ure 1 so that photomicrographs could be obtained aver a vertical trans- 
verse of the spray cross section. 


Drop-size and drop-velocity data were obtained by photographing the 
spray at distances of 1, 52; 14, and 18 inches downstream of the fuel 
injector center line. А camera unit developed at the NACA Lewis lab- 
oratory (ref. 8) was positioned at the 2-inch-square plastic windows, 
as shown in figure 1, to obtain photomicrographs of the isooctane spray 
in natural suspension in the high-velocity air stream. The object plane 
was centered along the test section center line (position 4, fig. 1) and 
also moved horizontally to obtain a horizontal traverse of the spray 
cross section (positions 3 and 5, fig. 1). Thus, by adjusting the cam- 
era unit in a horizontal direction and the fuel injector in a vertical 
direction, a seven-position traverse was made of the spray cross section 
at each of the four distances downstream of the injector. 


The lens system used in the camera unit had a magnification factor 
of 21, a l-centimeter-square field area, a field depth of 1 millimeter, 
and resolved droplets as small as 5 microns in diameter. Stopped images 
of the high-velocity droplets were obtained on the film, as shown in 
figure 2, by means of the rotating mirror and synchronized flash system 


shown in figure l. 
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The original droplet camera described in reference 1 was modified 
considerably in this investigation, since the additional factors of high 
droplet concentrations and varying droplet velocity were considered. 
The three-phase 400-cycle inverter was rewired and a rheostat was con- 
nected in the field circuit so that the speed of the rotating mirror 
could be varied from 3,000 to 11,000 rpm. The pair of 0.05-microfarad 
condensers were replaced with a pair of 1.0-microfarad condensers which 
gave & higher intensity flash when discharged across the two 0.014- by 
O.l-inch magnesium-ribbon electrodes spaced 3/8 inch apart and sand- 
wiched between split cylinders of bakelite insulation. Also, the lens 
pystem was redesigned and a modified aerial camera with the shutter and 
lens removed was used primarily as a holder for the 9.5-inch-width roll 
film. 


By using 1.0-microfarad condensers, which gave a very brilliant 
flash, it was possible to test various fuel dyes, light filters, film 
types, and developing solutions for obtaining high-contrast negatives. 
Best results were obtained with a clear fuel (not dyed), a Wratten red 
filter No. 25 mounted on the light source, and a shellburst panchromatic- 
type film developed with & D-11 solution. The red light, because of its 
relatively long wave length, gave & minimum of light scattering in the 
Spray, апа thus sharp droplet images were obtained on the high-contrast 
shellburst film. The D-11 solution was a standard developer generally 
used for high-contrast work. 


The following procedure was used in obtaining photomicrographs of 
the spray: The air velocity in.the test section was set at 140 or 180 
feet per second. An air temperature of 829239 F and static pressure of 
29.520.353 inches of mercury were used in both cases.  Isooctane metered 
at a flow rate of 100 pounds per hour with a rotameter was injected from 
the contrastream orifice into the air stream with a pressure drop of 55 
pounds per square inch across the orifice. Photomicrographs of the 
spray were then taken over a range of mirror speeds. From the photo- 
micrographs, the correct mirror speed required for obtaining stopped- 
droplet images on the film and also droplet velocity was determined. 
Photomicrographs of the spray were taken at the proper mirror speed and 
analyzed to obtain drop-size-distribution data. An average of seven pic- 
tures were obtained at each of the seven points in the traverse of the 
spray cross section, and all the droplets in focus and stopped on the 
film were measured to obtain at least 1400 droplets for each analysis. 
The same procedure was repeated at distances of 1, 55, 14, апа 18 inches 
downstream of the injector to obtain vaporization-rate data. 


Measurements of droplet diameters were made from the photomicro- 
graphs with an eyepiece having small-scale subdivisions of 0.10 milli- 
meter which gave droplet diameters accurate within +3 microns. The 
camera was calibrated with a Bausch and Lomb microscope slide having 
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small-scale subdivisions of 10 microns. Since 1400 to 1800 drops were 
measured for each analysis of the spray cross section, a statistical 


error factor of approximately 3= percent was calculated for each cummu- 


lative distribution analysis. 


The wet-bulb temperature for isooctane droplets evaporating under 
the conditions of this investigation was obtained experimentally. From 
the data on drop diameter, velocity, and wet-bulb temperature, it was 
possible to calculate vaporization rates and drag coefficients for the 


isooctane sprays. 


RESULTS AND ANALYSIS 
The following steps were employed in the analysis of the data: 


(1) A spray-vaporization-rate equation based on the mean drop 
diameter Deo was written for the case of simultaneous acceleration 


and evaporation of sprays. 


(2) Drop-size-distribution data, obtained from photomicrographs, 
were used to calculate the mean drop size 120 for each spray analysis. 


(3) Drop velocities were calculated and used to determine vaporiza- 
tion times and droplet drag coefficients. | 


(4) A comparison was made between the predicted and observed max- 
imum droplet diameter in the spray. 


(5) Spray vaporization rates based on mean drop diameter 7020 
were compared with values predicted by single-drop vaporization-rate 


equations (refs. 5 and 6). 
(6) The weight percent of the spray evaporated at distances of 1, 


55, 14, and 18 inches downstream of the injector was determined and 
compared with values predicted by an empirical expression (ref. 4). 


Vaporization-Rate Equations 


The instantaneous vaporization rate dm/d@ of a drop may be de- 
termined from the heat-transfer equation 


dm 
ag Њу = ВАДЕ = xDk,AtNu . (1) 


` 


СА-2 
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when the change in sensible heat of the liquid is negligible compared 
with latent heat requirements. In recent investigations (refs. 5 and 
6), the following expression was obtained for the heat-transfer Nusselt 
number: 


0.6 0.5 
к – DD e / pAuD gl =) 
Nu = — = 2 + 2.58XL0OY | ——— 2 
Kg | юс) Ky (2) 


where the product of the Reynolds and Schmidt numbers is given as a 
single dimensionless number equal to одар Љ w? and the momentum- 


transfer ratio gi /c* is the ratio of gravitational to viscous forces. 


Equation (1) may be rewritten by expressing the instantaneous 
vaporization rate in terms of the change in drop агеа with time and 
solving for the Nusselt number as follows: 


=: PTH, ар? (5) 


M = V 
Akt аө 


where Nu is the heat-transfer Nusselt number for a single drop. 


The change in area of the spray A may be written as the sum of 
(1) the area of the completely evaporated drops, and (2) the change in 
area of the partially evaporated drops: 
D= JAAGNu Ређа 
AA = x > | пре + xAAO 2-2 n Nu (4) 


D= „ЈАДОМа. 


where А = 4k At/prH,- During each experimental time increment ЛӨ, the 


decrease in spray area due to complete evaporation of small drops was 
negligible (approximately 1 percent) compared with the decrease in area 
due to partial evaporation of the larger drops in the sprays. Thus, 
the &bove expression was rewritten as 


D=D 
АА a Dn Na (5) 
D=0 


For this investigation, equation (2) may be written 
Nu = 2 + 0.39 ReÜ-8 


(since Ъ "E 2.7x107* g/(cm) (sec), gi /c* = 1.024x1071¢, апа 
(к/к) 0:2 = 1.64.) By combining equations (4) and (5) and dividing 


D= 
each side of the resulting equation by тор п, the expression 


e p GS ра. ae eS а — —— — - КЕЧ =: 
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0.6 


2 

AD Aup 
2-5 М2 + 0.39 = 
ЛӨ На 


D=Dn Dn 
D>, °F D а (6) 
D=0 D=0 


is obtained, where D8 = Хард ха. It was found experimentally that: 


= 0.90 + 0.02 (7) 


0-6 
D26 


Thus, combining this expression with equation (6) and passing to the 


limit yields 
aD So aoe, \ "P kt 
айо AD + oss (ан) ји (e 


where the term in brackets represents the Nusselt number based on the 
mean drop diameter Dog. 


The mean Nusselt number №20 may be determined from equation (8) 
апа compared with the single-drop Nusselt number given by equation (2) 
when the spray parameters авдбо, Au, and “At are known. Since the 


surface temperature of each drop could not be determined experimentally, 
the wet-bulb temperature of isooctane was obtained experimentally and 
assumed equal to the surface temperature of a single drop having a mean 
drop diameter equal to 020. The remaining parameters 40550 and Au 


were determined experimentally. 


Drop-Size-Distribution Analyses 


Determination of mean drop sizes. - Analyses of the drop-size dis- 
tribution for isooctane sprays at distances of l, 52, 14, and 18 inches 


downstream of the fuel injection station were made for air-stream veloc- 
ities of 140 and 180 feet per second. Table I shows the drop-size count 
obtained at each of the seven positions in the spray cross-section tra- 
verse for & representative test. А complete traverse of the positions 
1 to 7 and total drop-size count are included in this table. Eight 
total-drop-size counts are recorded in table II, and a typical plot of 
drop-size distribution is shown in figure 5. 


~ — ——— 08 не јаја чаран чи a ои 
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Three of the most common mathematical expressions for drop-size 
distribution are the Rosin-Rammler, Nukiyama-Tanasawa, and log- 
probability equations. A description of the evaluation and use of these 
expressions is given in reference 9. They may be written as follows: 


Rosin-Rammler , 


GR _ -99-1 ,-(D/5)* (9) 
dD үр 
Nukiyama-Tanasawa, 
6 
dR р 5 -bD 
47120296 | (10) 
Log-probability, 
2 
dv = -5. e 9 у? (11) 


бу Ws 


Drop-size-distribution data were plotted for each of the eight 
analyses using each of the three distribution equations. Samples of 
these plots are shown in figures 4 to 6. The Nukiyama-Tanasawa and 
log-probability expressions gave good agreement with values of Рго ob- 
tained by direct integration of the data. Results for the test condi- 
tion Us = 140 feet per second and x = 5.5 inches are shown in the 
following table: 


Mathematical expression Mean drop diameter, 
D202 
microns 


Rosin-Rammler 
Nukiyama-Tanasawa 
Log-probability 
Direct integration 


The Rosin-Rammler expression gave poor agreement since q in equation 
(9) was found to be approximately equal to 3. When q is equal to 3 

in this expression, number mean diameters such as D20 are predicted 
to be equal to zero. This, of course, limits the usefulness of equation 
(9). A similar limitation of the Rosin-Rammler equation was found in 


reference 9. 


Comparison with predicted mean drop sizes. - Reference 10 gives 
the following empirical expression for predicting the Sauter mean diam- 
eter Dzo produced by а gas atomizing nozzle: 


— — — M А —À — 75 7-- 


12 NACA TN 3265 


0.45 1:9 
Пар = 1920, r3 + 597 [12 (2225 (12 
“бв er сог. Qs 


For this investigation, 9/9, is small, and the expression reduces to 


= 2220 "т = 221218 (15) 


for isooctane. Ву using the Nukiyama-Tanasawa size-distribution ехргев- 
sion, Пао may be related to Роб as follows (ref. 11): 


Dao = 0.695 Dao 


which may be substituted into equation (13) to give 


7288 
D20 = зүс (14) 
8 


Values for the mean drop size Dog obtained l inch downstream of 


the injector for the two air-stream velocity conditions are listed in 
the following table, with values obtained from equation (14) for 
comparison: 


Air-strean Mean drop size, ор» microns 
velocity, 


Us, Experimental From eq. (14) 
ft/sec 


The agreement is good at the higher air-stream velocity. However, the 
effect of U, on Dog appears to be greater for the calculated values 


shown in the preceding table. This may be because of the different 
type of nozzle used in the Nukiyama-Tanasawa atomization studies. It 
is, therefore, doubtful that the Nukiyama-Tanasawa equations can be 
used to predict initial drop size over wide ranges of operating condi- 
tions with the simple orifice injector used in this investigation. 


Analysis of Drop-Velocity Data 


Calculation of drop velocity. - The velocity of droplets which ap- 
peared on the photomicrographs as stopped images was calculated from the 
expression 


4712) 
Ug = М 


(15) 
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Sharply defined droplet images were obtained on the film when mirror- 
speed settings gave values of drop velocities calculated by the preced- 
ing expression within +5 feet per second of the actual drop velocity. 
Since the duration of the magnesium flash was 4 microseconds, a droplet 
travel of 6 microns was permissible and droplet images on the film 
could be measured within +5 microns. Droplet images were badly smeared 
and lacked sufficient contrast for measuring when mirror-speed settings 
gave drop velocities which were not within +5 feet per second of the 
drop velocity. 


The proper mirror speed was determined by taking photographs at a 
low mirror speed and gradually increasing the speed until stopped drop- 
let images were obtained. The photographs showed that droplets of all 
sizes (table I) were moving at approximately the same velocity, because 
droplets of all sizes were stopped at a particular value of mirror speed. 
It was also established experimentally that a range of mirror speeds 
could be used in which photomicrographs were obtained showing stopped 
images of drops of all sizes at each of the mirror speeds. The range 
of mirror speeds resulting in definition of the drops was greatest at 
the station 1 inch downstream of the fuel injector; at this location the 
range of mirror speeds giving drop definition was about 20 percent of 
the mean value. Approximately the same drop-size distribution and mean 
drop Size were obtained from the highest and the lowest mirror speeds 
giving drop definition. This range of mirror speeds was probably due 
to the liquid jets penetrating several inches into the air stream and 
being broken up into drops in the penetration zone shown in figure 7. 
Thus, droplets formed at the maximum point of penetration had time to 
accelerate to & relatively high velocity, compared with drops formed at 
the orifice, before reaching the camera station. Images of droplets of 
all sizes could be stopped at a particular value of mirror speed with 
the injector at each of the different locations. 


Minimum mirror speeds which gave stopped droplet images were used 
in obtaining photomicrographs for each spray analysis so that the in- 
jector orifice could be considered to be the point of initial formation 
and acceleration of the spray. Drop velocities are plotted against the 
distance downstream from the fuel-injector orifice in figure 8 for the 
droplets starting at zero initial velocity at the injector orifice. For 
the case of solids injected into the 180-foot-per-second air stream, a 
5-micron-diameter sphere would accelerate to approximately stream veloc- 
ity while traveling a distance of 1 inch, Whereas 5-micron drops were 
found to have a minimum velocity of 85 feet per second, as shown in fig- 
ure 8. Thus, the acceleration of small drops was found to be consider- 
ably lower than that predicted for small solid spheres. 


Calculation of vaporization time. - Vaporization time was са1си- 
lated by first determining an empirical relation between drop accelera- 
tion and velocity difference, and then integrating the expression in two 
successive steps to obtain an equation giving distance of droplet travel 


— ee ee 2 + m 
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as a function of stream velocity and time. Acceleration was determined 
from a plot of drop velocity squared against distance (fig. 9) and the 
acceleration equation 


2 2 
„Ча gg. 
~ Bx 


Thus, the acceleration at any distance x was equal to one-half the 
slope of the curve in figure 9. 


In order to determine a relation between drop acceleration and 
velocity difference, the plot shown in figure 10 was prepared. From 
this plot, the following empirical expression was obtained: 


а = 333 (Ug - Ug) (16) 


Since acceleration may be expressed as the change in drop velocity with 
time 04/20, the preceding expression may be written ав 


Ua Ө 
= 355 | аө 


- Оа О 


Integrating between the limits shown gives 


U =U, (1 - е-5550) (17) 


By substituting 4х/40 for droplet velocity in the preceding ex- 


pression, 
x Ө Ө 
dx = Ug 1 49 - U, е-5350 ag 
O О О 


Integrating gives the final expression 


6. 
1-е 
а lo - 2-7 (18) 


is obtained. 


where the vaporization time Ө is an implicit function of the stream 
velocity U, and the distance x. 


ee эи een ДИ ЧЕҸЕ 
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Calculation of drop drag coefficients. - The following force bai- 
ance for an evaporating droplet accelerating in an air stream was de- 
rived (appendix A): 


Е = ша + Au дш/49 (19) 


The instantaneous aerodynamic pressure force of the air stream accelerat- 
ing the vaporizing drop may be expressed as 


F = = pgAt (Ли) “Ср (20) 


The following equation for the drop drag coefficient was obtained 
by combining equations (19) and (20) and substituting into the expres- 
sion terms for drop volume, drop frontal area, and vaporization rate 
(given by eq. (12)): 


8k AtNu 
cp = 5 -— + “Бү (21) 
Ре (Ли) vPg 


(For this investigation, py, = 0.69 g/cu сш, Pg = 12x1074 g/cu cm, 
Kg = 5.9410-5 &-са1/ (ст) (вес) (0С), At = 259 C, and H, = 75 g-cal/g.) 


Drop drag coefficients were calculated from the preceding expression 
and plotted against the drop Reynolds number as shown in figure 11 for 
five velocity differences. This figure shows that the droplet drag 
coefficient for evaporating liquids cannot be correlated directly with 
the Reynolds number to obtain a single curve for all velocity differ- 
ences.  Drag-coefficient values for solid spheres in air streams (ref. 
12) are also plotted in figure 11. 


An empirical expression was obtained for the drop drag coefficient 
by combining equations (16) and (21) to give 


Oy = 444 P D Вк АСНЫ 


е (22) 
p, Au H P PAu 

where the constant 444 contains the units reciprocal seconds required 
to balance the equation dimensionally. Further tests on the accelera- 
tion of drops are needed using other liquids and varying air-stream 
pressure and temperature conditions in order to establish a more general 
expression. 


Comparison of predicted and observed maximum droplet diameters. - 
When the aerodynamic forces over an accelerating drop exceed the surface 
tension forces, the drop will shatter. Thus, it is possible to calcu- 
late the maximum drop diameter that will be found in a spray under the 
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condition of maximum acceleration at the injector orifice. Some investi- 
gators (refs. 13 and 14) have studied the break-up of water drops, and 
equated aerodynamic pressure forces and surface tension forces as 


follows: 
l 2 
Cp = Pg(Au)“A' = ло 


where А! and D are the maximum drop frontal area and diameter, re- 
spectively. Maximum drop diameters, calcul&ted from this expression, 
are compared with observed values in the following table, which also 
shows values of the critical Weber number (We. = раџер /20): 


Air-stream | Maximum drop diameter, | Critical Weber 


velocity, number, Weges 
бы, based on ob- 
ft/sec served D, 


140 127 126 
180 108 112 


Thus, agreement was obtained between values of Dy calculated for at- 
omization at the orifice and observed values obtained 1 inch downstream 
of the orifice. 


Fuel-Spray-Evaporation Calculations 


Comparison of fuel-spray and single-droplet vaporization rates. - 
From the drop-size distribution, drop velocity, and wet-bulb data, the 
instantaneous vaporization rate of the spray expressed as the change in 
area of the mean drop size Doo with time may be obtained by plotting 


ро, against 0, as shown in figure 12. In order to compare single- 


droplet vaporization-rate equations with experimental spray vaporization 
rates, equations (2) and (3) were used to calculate 4054/40 fora 


single droplet, and the results were plotted in figure 12. Experimental 
and predicted values of Dog agree within approximately +5 percent. 


Thus, the assumption that the mean surface temperature is closely ap- 
proximated by the wet-bulb temperature appears valid in this investiga- 
tion. This validity was possibly a result of low over-all fuel-vapor 
concentration. Although the concentration of drops was high near the 
injector orifice, the velocity difference Au was also very high 
(droplet film thicknesses were small) so the spray evaporated essen- 
tially as a cloud of isolated droplets. Simplified vaporization-rate 
equations are derived in appendix B. 
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Comparison of predicted and experimental weight percent of spray 
evaporated. - Since the change in spray mass due to the complete evap- 
oration of small drops was negligible compared with the change in mass 
due to the partial evaporation of larger drops in the spray, the weight 
percent of the spray evaporated E was calculated from the expression 


(095); - (084), 


0%), 


Е = 100 


Values of 925 Obtained by direct integration of the drop-size- 


distribution data, are given in table III with the weight percent of 
the spray evaporated. Initial values of Dog were obtained from the 


curve in figure 12 at Ө = О. Using these values, initial values of 
рхо меге determined from the expression 


Рао = 1.10 Dog (23) 
Experimental data showed the ratio Dzo/D20 was equal to 1.1040.02, and 
the Nukiyama-Tanasawa expression for mean drop sizes gives 

Dzo = 1.13 D20 | 


Table ТТТ also gives values for the weight percent of spray evap- 
orated as calculated from the following empirical expression derived in 
reference 4: 


E | Ts)" Us ) 1.2 0.48 0.84 
100 -E = 9:5° 10 E Ps, Pp х (24) 


The agreement is good except for the station 1 inch downstream of the 
fuel injector, which represents an extrapolation of equation (24), since 
the sampling probe technique could not be used for distances of less 
than 5 inches from the fuel injector. 


SUMMARY OF RESULTS 


A droplet camera was used to obtain drop-size distribution and 
drop-velocity data for isooctane sprays. An analysis of heat-transfer 
and drag data for the sprays evaporating and accelerating in turbulent 
air streams gave the following results: 


1. Mean drop sizes were determined by direct integration of the 
data and compared with values obtained from an analysis of the data using 
the Rosin-Rammler, Nukiyama-Tanasawa, and log-probability expressions 
for drop-size distribution. Agreement was obtained with the Nukiyama- 
Tanasawa and log-probability expressions. 


- Б Бин м л dc a лз АЕ 5E 
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2. The mean drop diameter Dog predicted by the empirical Nukiyame- 


Tanasawa expression for a different type of nozzle agreed with the ex- 
perimental value of Dep obtained at a distance of 1 inch downstream 


of the injector in a 180-foot-per-second velocity air stream. The effect 
of stream velocity on drop size predicted in the Nukiyama-Tanasawa ex- 
pression is different, however, from that obtained experimentally. 


5. An empirical expression was derived for the drag coefficients of 
isooctane drops accelerating and evaporating in turbulent air streams. 


4. An expression, obtained by equating aerodynamic pressure forces 
and surface-tension forces for drops, was used to calculate the maximm 
drop size existing in the spray for the condition of maximum spray ас- 
celeration. Photomicrographs gave maximum drop diameters which were in 
agreement with this expression. | 


5. The vaporization rate of an isooctane spray based on the mean 
drop diameter Рос (defined as a drop having an area equal to the total 


area to total number ratio for all the drops in the spray) was found to 
correlate single-drop vaporization rates. The weight percent of the 
spray evaporated downstream-of the injector was determined and found 
to agree with values obtained using a different technique. 


CONCLUDING REMARKS 


Spray-vaporization~rate data obtained in this investigation corre- 
late single-droplet vaporization rates. The Reynolds number effect 
during spray acceleration was found to be very important in determining 
the time required to evaporate the spray. Approximately 50 percent by 
weight was evaporated during the acceleration period. Thus, when more 
data on spray acceleration and atomization are available, the ехргев- 
sions used in this investigation may be tested and extended in the field 
of spray atomization, acceleration, and evaporation. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 30, 1954 
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APPENDIX А 


CHANGE OF MOMENTUM OF EVAPORATING DROP 


А mass m + An is considered to be moving at a velocity Ug as 
shown for condition I: 


Condition I Condition II 
{ЧҮ > Ua > Џа 
Ча + АЏа 


The velocity of the mass m is increased to Ug + Ag, and the velocity 
of the element Ат is increased to stream velocity Us, as shown above 


for condition II. Thus, the change in momentum AG is given ав 
follows: 


AG = ш (Ug + Aa) + Am U, - (m+ Am) Ua 
which reduces to 
AG = mUa + Am(Us - Оа) 


Dividing by AO and passing to the limit yields the instantaneous force 
F as follows: 
aU 
d аш 
=m 45 + (Us - Ua) a 


ас 
r = 49 49 


where dm/d@ is the instantaneous vaporization rate. 


"— a чы «= тоте - Mm mortcm a PE чи e YS п A SE рна А a a Mh пати иг — санаан —- С - 
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FUEL-SPRAY-EVAPORATION CALCULATIONS 


The time required to evaporate a spray may be divided into two 
periods: (1) the period required to accelerate the spray to approxi- 
mately stream velocity, and (2) the remaining period of evaporation in 
Which the spray travels at approximately stream velocity. In this in- 
vestigation, the spray acceleration period was studied, and the spray 
vaporization rate (based on the mean drop diameter Dep) was found to 
agree with equation (8): 


2 
Би кл 
a + 0.55 ү 


“48 "р 


The final period of spray evaporation has been treated in reference 
7 for the spray traveling at approximately stream velocity. Thus, a 
complete analysis of fuel spray evaporation may be made by treating 
each period separately and combining the results. 


Evaporation During Spray Acceleration Period 


2 
Determination of change in Doo during spray acceleration period. - 
For this investigation, 


^g" per ме те iut; а db cx ава =“ = = 159х100 


t, Ор » . Ф » Ф » Ф . е Ф • * е е • • » » • LJ Ф Ф 25 


Es g/cu cm . » е е ч *. * е Ф e » € е Ф е о е . Ф е Ф Ф Ф Ф 0.69 
Bys g-cal/g • ө > е э o ө € . © • © + o o э e оо ө э o Өө э 15 
and the vaporization-rate equation may be written ав 
2 0.6 
Dao 20 
64 = 11.65x10-5| 2 ерес | —2 
аө Ha 


Equation (17) gives Ли in terms of 0, and Ө. By substituting this 
expression into the vaporization-rate equation, 


2 0.6 
dD 0 р 

a = 11.65х10-5|2 + еке e 200 9 (25) 
a 
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is obtained. This expression may be graphically integrated by taking 


small increments of Ө and calculating the change in 05. 


It was found experimentally that the sprays had accelerated to 
approximately stream velocity in 0.014 second. Although 20 decreased 
approximately 33 percent in this time, Dos. decreased only 5 percent. 
Thus, if Бо? is assumed approximately constant, the vaporization rate 
equation may be integrated to give 


| эд )49.0а| 0-6 
ApS, = 12.65x10"5{0 028 + al 2022 za | 


a 


= 0.014 ,-200 Ө о ET 67200 Ө - 
© 200 
0-0 


where С = 0.06 for the condition of complete spray acceleration. By 
using the following values: 


since, 


(025) 1? ва cm е * е . * • • • Ф » » » е » е е е е • • * e • 2100x1078 


Us; ft/sec € ө ьо © ө © Ф ө ө © ө өө өө a өө © > ө өр $4 Фө &© $9 @ 140 
На» poises e e 9" > ө ө ө ө e Ф > ө 9 өр > э ө © ө өө o ө où ө + 1 .8х10-& 
Da» gn сц CM . «© © © e * © © © . > o e э ө 9 © э ө э © ++ o ò ө 12 .04x1074 


the expression 


Ans, = (526 + 375)х10-8 = 701Х1078 ва cm 
may be written. Graphical integration of equation (25) gave р = 
100Х10-8 square centimeter (fig. 12) for О, = 140 feet per second, and 


Ө = 0.014 second. Thus, a simple expression was obtained, since the 
change in the Reynolds number was controlled primarily by the change 
in velocity difference Au and was affected very little by the change 
in 820. 


Calculation of weight percent evaporation during spray acceleration 
period. - The weight percent evaporated E during the spray accelera- 


tion period may be approximated by the expression 


E = 100 x р У СОДОН (26) 


since it was found experimentally that 050/020 = 1.1040.02. Values 


of E calculated from this expression and experimental values of E 
are compared in the following table: 
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Air-stream | Time, | Distance from 

velocity, injector, x, 
Us; in. Experimen- | Calculated 
ft/sec tal, from | from eq. (26) 


table III 


Evaporation During Period Following Complete Spray Acceleration 


The method of reference 7 may be used to calculate the weight 
percent evaporated after the spray has accelerated to stream velocity. 
In this method, the Rosin-Raumler expression (eq. (9)) is used. When the 
values of D and а are known, the weight percent unevaporated 100 - E 
may be determined. 


The relation between D and Dj шау be derived in the following 


manner: Тһе Nukiyama-Tanasawa equation for mean drop sizes (ref. 11) 
predicts Dzj/Dog = 1.29 which agreed within 2 percent of experimental 


values of аа /Оар obtained in this investigation. Also, it was shown 
in reference 7 that the Rosin-Rammler expression for mean drop sizes 


gave Р = м2, D3] ава close approximation for D when 2<n<4. Thus, 


| D = 1.82 0,0 
may be written. Figure 12 shows Doo = 43.6 microns at U, = 140 feet 


per second and x = 5.5 inches (@ = 0.0058 sec). By substituting this 
value of Dog into the preceding expression, D = 79.5 microns is ob- 


tained. This value agrees with figure 4 which shows D = 80 microns. 
Experimental values of D/Deg agreed with the constant 1.82 within 2 


percent. 


In this study, the constent n in equation (9) was approximately 
5 throughout the spray acceleration period. By using n = 5, and 
calculating D from the expression 


D U 0.6 0.5 
а - 11.65X10-? {0.028 + 1.175Х1079 [юне] | 


D = 1.82 (22), 


а, 


The weight percent evaporated Е following the acceleration period may 
be determined from the plot of reference 7 of 100 - E against „Ља 
(where Х = 23.3 1079. сте /вес and 9 is the time following the spray 
acceleration period). Thus, the total weight percent evaporated may be 
calculated as the percent evaporated during the spray acceleration period 
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plus the weight percent of the remaining fraction evaporated in the 
final evaporation period. 


Application of Method of Reference 7 to Spray Acceleration Period 


If it were assumed that the spray traveled at stream velocity 
(Au = О and Nu = 2) in the 140-foot-per-second air stream for the 
distance x = 18 inches (Ө = 0.014 sec), the weight percent evaporated 
would be calculated in the following manner: From the plot of refer- 
ence 7 of weight percent unevaporated spray (100 - E) against 4A0/D, 
the expression 


-5 
e/o 2 25.5X]10 "- X мон = 0.198 
85.5Х107 


may be obtained, from which 100 - Е = 88 percent, or E = 12 percent 
evaporated. This value of E is approximately 30 percent less than 
the experimental value given in table III. Thus, a very serious error 
would be made if the evaporation of the spray during the acceleration 
period was assumed the same as evaporation in a quiescent atmosphere. 
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TABLE Т. - DROP-SIZE-DISTRIBUTION DATA FOR SEVEN-POINT TRAVERSE 


QF SPRAY CROSS SECTION 


[Air-stream velocity, Од, 140 ft/ Sec; distance from injector 
orifice, x, 18 їп] 


Number of drops, п 


Traverse position? Total |Second |Total 


а tp 
"HERES 


2.5 ol o| i| 2 

7.0 2| 1| 5| 2 

12.0 6| 01:8| 7 

17.0 11 

21.5 21 

26.0 18 

51.0 6 9 

56.0 11 11| 15| 8 

40.5 10 4| 9| 8 

45.0 ii| 91-51 4| 2| 8 

50.0 8| 4| 4| 4| 5| 3| 7 

55.0 5| 6| 4| 3| 9| 1| 5 

59.5 2| 2| 11:41 4] 1| 4 

64.0 9| 2| 1| 3| 4| a] 5 6 

69.0 6l 4| 2| 3| 5| 011 6 

74.0 5| 2| of 2| ol ol o 10 

78.5 11 11 2| 3| 51 2! 1 О 

85.0 ї| ol 11.51 2| O| o 2. 

88.0 2| 2| 1| of 5! ој о 1 

95.0 01 2| oi o| 3] o| o 5 

97.5 ol ol of 31 ofl 1| o 1 
o| 1| ol Ol LI] 01 o О 
0| 1| ol ol о| о| o О 
ol 1| ol 01 Ol ol o О 


[ Tete раз |ізо iso |ізо so [iso iso] sos | 564 _ 


BAD = 4.75 microns. 
кз, 1% 


26 


NACA TN 3265 


Distance from injector, x, in. 


Number of drops, n 
8 12 10 10 16 20 25 26 


18 
11 
12 
8 
6 
4 
5 
3. 
1 


"л:озн атњо 


, | 
PY гоо о (л со (а 3 


FO FRFWUrF © 
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TABLE IIT. - WEIGHT PERCENT OF SPRAY EVAPORATED 


Distance D39X1010 P Weight percent evaporated 
from m Experimental | Calculated from 
injector , equation (23) 
х, 
in. | 


Air-stream velocity, U,, 140 ft/sec 


8calculated from Dzo/ D20 = 1.10 where Do) was 
obtained from fig. 12 at х = О. 


21 
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Figure 1. - Schematic diagram of test section and camera unit E 
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‘Figure 2. - Photomicrograph of isooctane spray in 180-foot-per-second air stream. 
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Figure 3. - Droplet-size distribution. 
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- Rosin.and Rammler analysis. 


velocity, 140 feet per second; distance from injector, 


5.5 inches; D, 80 microns; mean drop diameter Doo, 


21.1 microns; 


Figure 4. 
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Drop diameter, D, microns 


0.0837. 


-3 microns; b, 


Air-stream velocity, 140 feet per second; 


distance from injector, 5.5 inches; mean drop diameter Doo 41 


Figure 5. - Nuktyama and Tanasawa analysis. 
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Figure 6. ~ Log-probability analysis. Air-stream velocity, 140 feet per second; distance from 
injector, 5.5 inches; D, 66 microns; 8, 1.712; mean drop diameter Про; 40.7 microns. 
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Figure T. - Break-up of liquid fuel jet in turbulent air stream. 
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Figure 8. ~ Velocity of isooctane droplets at given distances dowetream of fuel-injector orifice. 
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Figure 9. - Acceleration of isooctane droplets. 
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Figure 10. - Relation between drop acceleration and 
velocity difference. 
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